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Abstract. We consider a six-parameter family of the square integrable wave functions for the 
simple harmonic oscillator, which cannot be obtained by the standard separation of variables. They 
are given by the action of the corresponding maximal kinematical invariance group on the standard 
solutions. In addition, the phase space oscillations of the electron position and linear momentum 
probability distributions are computer animated and some possible applications are briefly discussed. 
A visualization of the Heisenberg Uncertainty Principle is presented. 



The purpose of this Letter is to elaborate on a "missing" class of solutions to the time-dependent 
Schrodinger equation for the simple harmonic oscillator in one dimension. We also provide an 
interesting computer-animated feature of these solutions — the phase space oscillations of the 
electron density and the corresponding probability distribution of the particle linear momentum. 
As a result, a dynamic visualization of the fundamental Heisenberg Uncertainty Principle 27 is 
given 



45 , 62 



1. Symmetry and Hidden Solutions 

The time-dependent Schrodinger equation for the simple harmonic oscillator, 

2#, + 7/^,,. - x^V; = 0, (1.1) 
has the following six-parameter family of square integrable solutions 

i[a(t)x'^+5{t)x+K{t))+i(2n+l)-t{t) 

t) = , , ^ e-(^W-+^W) (/3 {t)x + e (t)) , (1.2) 

V2"n!/i (t) 

where (x) are the Hermite polynomials [sT] and 



fi{t) = /iQY/3gsin t -|- (2q;o sin t -|- cos t) , (1.3) 
«ocos2t + sin2t (/3^ + 4«2_iy4 

= M • 2^ , • ^ , ^ — ' ^^-^^ 

Po sm t + [2ao sm t + cos t) 

nt) = ^° (1.5) 

/3q sin^ t + (2q;o sin t + cos t) 

1 /3q sin t 

1 [t] = 7n arctan , (1.6) 

'° 2 2aosint + cost' ^ ^ 
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Sit) 
e{t) 

Kit) 



60 {2ao sint + cost) + eo/^o sint 

/Sq sin^ t + (2ao sin t + cos t)^ 
£0 (2ao sin t + cos t) — /3q5o sin t 



Pq sin^ t + (2ao sin t + cos t)^ 



Kq + sin^ t 
1 

+ - sin2t 
4 



/3q sin^ t + (2ao sin t + cos ty 

2 /o2 _ r2 



(1.7) 
(1.8) 

(1.9) 



/3q sin^ t + (2ao sin t + cos t) 



2 • 



(Here, /ig > 0, ao, /Sg 7^ 0, 79, ^o, ^^o, ^0 are real initial data.) These "missing" solutions can be 
derived analytically in a unified approach to generalized harmonic oscillators (see, for example, [o], 
[10] , [39] and the references therein). They are also verified by a direct substitution with the aid 
of Mathematica computer algebra system 33 , 45 , |62]. (The simplest special case /ig = /Sq = 1 
and ao = 'Jq = 6q = Sq = kq = reproduces the textbook solution obtained by the separation of 



variables 21 , 23 , 38 , 47 ; see also the original Schrodinger paper [55j; and the shape-preserving 
oscillator evolutions occur when ao = and (3q = 1. More details on the derivation of these formulas 
can be found in Refs. 44 and 50 ; see also the references therein.) 



On the other hand, the "dynamic harmonic oscillator states" (1.2)-(1.9) are eigenfunctions 

1 



E(t)^„ {x,t)-- 

of the time-dependent quadratic invariant, 

{p — 2ax — 6) 



n + 



E{t) 



+ (/3x + ef 



1.10) 



:i.ii) 



[a it) a) it) + (t) a (t)] 



dt 



(E) = 0, 



with the required operator identity [15], [54]: 

dE 



1 



+ t-'[E,H]=0, H=-{p'' + x^) 



dt ' " " 2 

Here, the time-dependent annihilation a (t) and creation a"'' (t) operators are explicitly given by 



1.12) 



a{t) 



1 
71 



f3x + e + i 



P 



2ax — 6 



W{t) 



1 
71 



I3x + e — i 



p — 2ax — 6 



;i.i3) 



with p = i ^d/dx in terms of our solutions (1.4)-(1.9). These operators satisfy the canonical 
commutation relation, 

a{t)a^t)-a^t)a{t) = l, (1.14) 



and the oscillator-type spectrum (1.10) of the dynamic invariant E can be obtained in a stan- 
dard way by using the Heisenberg-Weyl algebra of the rasing and lowering operators (a "second 
quantization" |l], [42], the Fock states): 

a{t)'^n{x,t) = ^^n-i{x,t), {t)^n{x,t) = VnTl^n+i{x,t). (1.15) 
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Here, 



(1.16) 

is the relation to the wave functions (1.2) with (t) = — {2n + 1) 7 {t) being the nontrivial Lewis 



phase 42 , 54 



This quadratic dynamic invariant and the corresponding creation and annihilation operators for 
the generalized harmonic oscillators have been introduced recently in Ref. |54j (see also 10 , 15 , 59 
and the references therein for important special cases). An application to the electromagnetic-field 
quantization and a generalization of the coherent states are discussed in Refs. l35] and 40 



The key ingredients, the maximum kinematical invariance groups of the free particle and harmonic 
oscillator, were introduced in |3|, |4), |25), |30|, |49| and |50| (see also (7|, |3l], |i8], |53 



57 , 58 , 63 



and the references therein). We establish a (hidden symmetry revealing) connection with certain 
Ermakov-type system which allows us to bypass a complexity of the traditional Lie algebra approach 

(A general procedure 



see 



17 



|41| and the references therein regarding the Ermakov equation) 
of obtaining new solutions by acting on any set of given ones by enveloping algebra of generators 
of the Heisenberg-Weyl group is described in [15].) In addition, the maximal invariance group of 



the generalized driven harmonic oscillators is shown to be isomorphic to the Schrodinger group of 



the free particle and the simple harmonic oscillator 44 , 49 , 50 



2. Discussion 



Quantum systems with quadratic Hamiltonians 
651 



20 , 26 , 46 , 64 



(see, for example, |2), |5), |6|, [To], (T4|, (Ts), (18 

[66j , [68j and the references therein) have attracted substantial attention 
over the years because of their great importance in many advanced quantum problems. Examples 
are coherent and squeezed states, uncertainty relations. Berry's phase, quantization of mechanical 
systems and Hamiltonian cosmology. More applications include, but are not limited to charged par- 
ticle traps and motion in uniform magnetic fields, molecular spectroscopy and polyatomic molecules 
in varying external fields, crystals through which an electron is passing and exciting the oscillator 
modes, and other mode interactions with external fields. Quadratic Hamiltonians have particular 
applications in quantum electrodynamics because the electromagnetic field can be represented as a 
set of generalized driven harmonic oscillators 



13 , 20 



The maximal kinematical invariance group of the simple harmonic oscillator 50 provides the 



six-parameter family of solutions, namely (1.2) and (1.3)- 



data (of the corresponding Ermakov-type system 17 , [39 



1.9), for an arbitrary choice of the initial 
41 , (44)). These "hidden parameters" 



usually disappear after evaluation of matrix elements and cannot be observed from the spectrum. 
How to distinguish between these "new dynamic" and the "standard static" harmonic oscillator 
states (and which of them is realized in a particular measurement) is thus a fundamental problem. 



At the same time, the probability density \'ilj{x,t)\ of the solution ( |l.2[) is obviously moving 
with time, somewhat contradicting to the standard textbooks 



21 , 23 , 38 , 47 , 55 , - an el- 



ementary Mathematica simulation reveals such space oscillations for the simplest "dynamic os- 
cillator states" 44 , |45] (see Appendix A for the Mathematica source code). 



|45| (see Appendix A for the Mathematica source code). The same is true 
for the probability distribution of the particle linear momentum due to the Heisenberg Uncer- 
tainty Principle 27 . These effects, quite possibly, can be observed experimentally, say in Bose 
condensates, if the nonlinearity of the Gross-Pitaevskii equation is turned off by the Feshbach reso- 
nance 11 , 19 , 32 , 52 , 56 , 60 . A more elementary example is an electron moving in a uniform 
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magnetic field. By slowly changing the magnetic field, say, from an initially occupied Landau level 
with the standard solution |38], |43], one may continuously follow the initial wave function evolution 
(with the quadratic invariant) until the magnetic field becomes a constant once again (a parametric 
excitation; see, for example, lo], 15 , 37 , 46 and the references therein). The terminal state 



will have, in general, the initial conditions that are required for the "dynamic harmonic states" 



(1.2)-(1.9) and the probability density should oscillate on the corresponding Landau level just as 



our solution predicts. However, it is still not clear how to observe this effect experimentally (but 
these "dynamic harmonic states" will have a nontrivial Berry's phase [5], [6], 54 



61 



One may imagine other possible applications, for example, in molecular spectroscopy 

[66j , and cavity quantum electrodynamics [12j, [13j 



of crystals, quantum optics 24 



22 



62]). 

46], theory 



35 , 67 



We believe in a dynamic character of the nature |28j. All of that puts the consideration of this 
Letter into a much broader mathematical and physical context — This may help better understand 
some intriguing features of quantum motion. (Our example shows that the separation of variables 
for the time-dependent Schrodinger equations may not always give us the "whole picture".) 

Acknowledgments. We wish to thank Professor Sir Michael Berry, Professor Andrew Bremner, 
Professor Carlos Castillo-Chavez, Professor Victor V. Dodonov, Professor Martin Engman, Pro- 
fessor Georgy Th. Guria, Dr. Christoph Koutschan, Dr. Sergey I. Kryuchkov, Professor Elliott 
Lieb, Dr. Francisco F. Lopez-Ruiz, Professor Alex Mahalov, Professor Vladimir I. Man'ko, Dr. 
Benjamin R. Morin, Professor Peter Paule, Dr. Audrey M. Puchkov, Priv.-Doz. Dr. Andreas 
Ruffing, Professor Simon Ruijsenaars, Professor Vladimir M. Shabaev, Professor Erwin Suazo, Dr. 
Nikolay Tishchenko, Professor Luc Vinet and Professor Doron Zeilberger for support, valuable 
discussions and encouragement. This paper has been initiated during a short visit of one of the 
authors (SKS) to The Erwin Schrodinger International Institute for Mathematical Physics and we 
thank Professor Christian Krattenthaler, Fakultat fiir Mathematik, Universitat Wien, for his hos- 
pitality. This research is supported in part by the National Science Foundation-Enhancing the 
Mathematical Sciences Workforce in the 21st Century (EMSW21), award # 0838705; the Alfred P. 
Sloan Foundation-Sloan National Pipeline Program in the Mathematical and Statistical Sciences, 
award # LTR 05/19/09; and the National Security Agency-Mathematical & Theoretical Biology 
Institute — Research program for Undergraduates; award # H98230-09-1-0104. 



Appendix A. Mathematica source CODE lines 45 



Example 1 



The following animation is for the dynamic ground state n = 0, using ao = 7o = = 0, /3o = 2/3, 
<5o = 1: 



in[i]:= Animate 



Plot 



72(x- Sin[^-.(-l+T)]y 

9\/2 e cos[3io"(-i+T)] 



97 + 65 Cos [^7r(-l + T)] 



, {a;, —3.5, 3.5}, AxesLabel -> 



{x, (Abs[V'])^2}, PlotRange {0, 2.3}, Filling ^ {1 ^ Bottom}, PlotStyle^ {Thick, Blue} 



{T, 1001} 
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Out[l]= 




(A) (B) (C) (D) (E) 

Figure 1. Subfigures (a)-(e) are a few stills taken from the mathematica movie an- 



imation 



45 . Starting with (a) they denote the oscillating electron density (blue) of 



the ground "dynamic harmonic state" of the time-dependent Schrodinger equation 



(1.1). These space oscillations complement (with the help of Mathematica) the cor- 
responding "static" textbook solution (clear) [55]. (The color version of this figure is 
available only in the electronic edition.) 

Example 2 



The fohowing animation is for the first excited dynamic state n = 1, using ckq = 7o 
/3o = 2/3, 6o = 1: 



0, 







Plot 


IK 



In [2]:= Animate 

/(97 + 65 Cos[^ni-l + T)]) 



72(x- Sin[^,r(-1 + T)])- 



500 



7r(-l +T) 



3/2 8 



e 9 X 



{x, -4.5, 4.5}, AxesLabel ->{x, (Abs[V'])^2}, 



Out[2] = 



PlotRange {0, 1.67}, Filhng ^ {1 ^ Bottom}, PlotStyle {Thick, Blue} 



{T, 1001} 





(A) (B) (C) (D) (E) 

Figure 2. Subfigures (a)-(e) are a few stills taken from the mathematica movie 



animation 45 . Starting with (a) they denote the oscillating electron density (blue) 
of the first excited "dynamic harmonic state" of the time-dependent Schrodinger 



equation (1.1). These space oscillations complement (with the help of Mathematica) 
the corresponding "static" textbook solution (clear) [55^. (The color version of this 
figure is available only in the electronic edition.) 
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Example 3 

The following animations simultaneously show the phase space oscillations of the electron density 
and the momentum probability distribution, according to the Heisenberg Uncertainty Principle, for 
the dynamic ground state n = with parameters uq = 'Jq = Eq = kq = 0, I3q = 2/3 and 5o = 3/2: 



In [3]:= Animate 



Plot 



ls(^-2x+3 Cos 




-l+T) 


r 


9V2e -™Cos 




-1+T) 





18(2a! + 3 Cos[5lg„(249+ T)])^ 
9\/2 e Cos[^^i-l + T)] 

97 + 65 Cos [^7r(-l + T)] 



., {j;, -4.5,4.5},AxesLabel-> {{x,p},{{Ahs['ilj])^2, 



^97-65 Cos[^ni-l+T)] 

(Abs[a])^2}},PlotRange {0, 2.3}, Filling ^ {1 ^ Bottom}, PlotStyle {Thick, Blue} 



{T, 1001} 



Out [3]= 




(F) 



(G) 



(H) 



(I) 



(J) 



Figure 3. Subfigures (a)-(j) are a few stills taken from the mathematica movie an- 
imation [45] . Starting with (a) they denote simultaneous oscillations of electron den- 
sity (blue) and probability distribution of momentum (clear) for the ground "dynamic 



harmonic state" of the time-dependent Schrodinger equation (1.1 ). These phase space 
oscillations complement (with the help of Mathematica) the corresponding "static" 
textbook solutions [23], [55]. (The color version of this figure is available only in the 
electronic edition.) 



One immediately recognizes from these animations that the particle is the most localized at 
the turning points when its linear momentum is the least precisely determined, as required by the 
fundamental Heisenberg Uncertainty Principle 27 — The more precisely the position is determined. 



the less precisely the momentum is known in this instant, and vice versa (see also [§]). In the creator 
own words — "If the classical motion of the system is periodic, it may happen that the size of the 
wave packet at first undergoes only periodic changes" (see Ref. 27 , p. 38). 
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According to (1.13), the corresponding expectation values are given by 

d 



(x) 



— [(2ao6o - Z^o'^o) sin t + eo cost] , 

d 



- [(2ao£^o - /3o'^o) cos t - Eq sin t] , 

Po 



dt 



{x) = {p), 



(A.l) 
(A.2) 



— (2q;o£:o ~ /^o^o) / f^o- This provides a classical 



with the initial data (a;)|^^Q = — ^o/Z^o {p)\t=o 
interpretation of our "hidden" parameters. 

The expectation values (x) and (p) satisfy the classical equation for harmonic motion, y" + y = 0, 
with the total energy 



{2aoeo - f3QSo f + el ^ 1 
2/3^ " 2 



For the standard deviations. 



{(Apy) = {p^)-{pf = {n+l) 



1 + 4ag + /3o + (4ag + (3^ - l) cos 2t - 4ao sin 2t 



((Ax)') = (x^) - {xf = {n+\) 



l+4a§ + /3o- [Aal + 131-1) cos 2t + 4ao sin 2t 



2/3g 



(A.3) 

(A.4) 
(A.5) 



one gets 



{{Apf){{/\xf) 



n + 



4/3 



(l + 4a2 + /34)2_ (|(|4^2^^4_^^^Qg2t-4aosin2t)' . (A.6) 



In the case of the Schrodinger solution 55 , when = 5q = Eq = Q and Pq = 1, we arrive at 
(x) = (p) = and 



{(Apf) = ((Ax)') = n + ^ 



(A.7) 



as presented in the textbooks 21 , |23], |24], [29], [38], [47]. The dependence on the quantum number 
n, which disappears from the Ehrenfest theorem 16 , [27j, is coming back at the level of the higher 
moments of the distribution. 



According to (A.6), 



((Ap)')((Ax)')= n + 



1 - 4al sin' 2t 



(A.8) 



provided that Aa^ + /Sq = 1; and the product is equal to 1/4, if n = and sin' 2t = 1. These are 
conditions for the minimum-uncertainty squeezed states of the simple harmonic oscillator (see, for 
example, 29 , [36]). For the coherent states ao = and /3q = 1, which describes a two-parameter 
family with the initial data {x)\^^q = —Eq and (p)lt=o ~ '^o- 

The corresponding wave functions in the momentum representation are derived by the (inverse) 
Fourier transform of our solutions (1.2) and (1.3)-(1.9) in AppendixB. Moreover, in Appendix C, we 
explicitly present the action of the Schrodinger group on the wave functions of harmonic oscillators 
and elaborate on the corresponding eigenfunction expansion for the sake of 'completeness'. 



More examples are available in the authors' websites. 
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Appendix B. The Momentum Representation 



For the wave functions in tfie momentum representation 

1 



a„ {p, t) 



(B.l) 



tlie integral evaluation is similar to Ref. [39j. As a result, the functions a„ (p, are of the same 

/So 



form (1.2)-(1.9), if ip^ — >• a„ and x — )■ p, with the initial data 



4«o + /3o 



2 ' 



/3i 



1 / 



5i 



2ao(5o + /3|^£o 



^1 



2q;o£o - /gp^o 

4ag + /3o 



(B.2) 

(B.3) 
(B.4) 

(B.5) 



The calculation details are left to the reader |62] (see, for example, Ref. 23 for the classical case) 



Appendix C. The Schrodinger Group for Simple Harmonic Oscillators 



The following substitution 



^i[a{t)x'^+5(t)x+K(t)) 



(C.l) 



where relations (1.3)-(1.9) hold, transforms the time-dependent Schrodinger equation (|l.l|) into 
itself with respect to the new variables ^ = P{t)x + e{t) and r = —7 {t) 



references therein). A Mathematica verification can be found in Refs. [33 



50 (see also 44 



and 62 . 



and the 



The eigenfunction expansion of the "dynamic harmonic states" with respect to the standard 



"static" ones can be obtain in an obvious way (see, for example, 37 and 43 for similar integral 
evaluations, the details will appear elsewhere). The corresponding matrix elements define the 
representation of the Schrodinger group acting on the oscillator wave functions. (The structure of 
the Schrodinger group in two-dimensional space-time as a semidirect product of SL (2, M) and Weyl 
W (1) groups is discussed, for example, in Refs. [T], [31] and [48].) 

An explicit time evolution of the (bosonic field) creation and annihilation operators for the "dy- 
namic harmonic (Fock) states" (with the embedded hidden Schrodinger group symmetry) can be 
easily derived from (1.13), (1.15) and (1.16). Applications to the quantization of electromagnetic 



fields are discussed in 35 . 



References 

[1] A. Akhiezer and V. B. Berestetskii, Quantum Electrodynamics, Interscience Publishers, New York, 1965. 
[2] V. Aldaya, F. Cossio, J. Guerrero and F. F. Lopez-Ruiz, The quantum Arnold transformation, J. Phys. A: Math. 
Theor. 44 (2011), 065302 (19 pages). 



HIDDEN SYMMETRY OF HARMONIC OSCILLATORS 



9 



[3] R. L. Anderson, S. Kumei and C. E. Wulfman, Invariants of the equations of wave mechanics. I, Rev. Mex. Fis. 

21 (1972), 1-33. 

[4] R. L. Anderson, S. Kumei and C. E. Wulfman, Invariants of the equations of wave mechanics. II One-particle 

Schrddinger equations, Rev. Mex. Fi's. 21 (1972), 35-57. 
[5] M. V. Berry, Quantal phase factors accompanying adiabatic changes, Proc. Roy. Soc. London, A392 (1984) 

# 1802, 45-57. 

[6] M. V. Berry, Classical adiabatic angles and quantum adiabatic phase, J. Phys. A: Math. Gen. 18 (1985) # 1, 
15-27. 

[7] C. P. Boyer, R. T. Sharp and P. Winternitz, Symmetry breaking interactions for the time dependent Schrddinger 

equation, J. Math. Phys. 17 (1976) #8, 1439-1451. 
[8] D. Cassidy, Uncertainty: The Life and Science of Werner Heisenberg , New York: W. H. Freeman, 1992; see also 

http : / / www. aip . org/history / heisenberg/ p08 . htm. 
[9] R. Cordero-Soto, R. M. Lopez, E. Suazo and S. K. Suslov, Propagator of a charged particle with a spin in 
uniform, m,agnetic and perpendicular electric fields, Lett. Math. Phys. 84 (2008) #2-3, 159-178. 
[10] R. Cordero-Soto, E. Suazo and S. K. Suslov, Quantum integrals of motion for variable quadratic Hamiltonians, 

Ann. Phys. 325 (2010) #9, 1884-1912. 
[11] E. A. Cornell and C. E. Wieman, Nobel Lecture: Bose-Einstein condensation in a dilute gase, the first 70 years 

and some recent experiments. Rev. Mod. Phys. 74 (2002) #3, 875-893. 
[12] V. V. Dodonov, Current status of dynamical Casimir effect, Physica Scripta 82 (2010) #3, 038105 (10 pp). 
[13] V. V. Dodonov, A. B. Klimov and D. E. Nikonov, Quantum phenomena in nonstationary media, Phys. Rev. A. 
47 (1993) # 5, 4422-4429. 

[14] V. V. Dodonov, L A. Malkin and V. L Man'ko, Integrals of motion, Green functions, and coherent states of 

dynamical systems, Int. .1. Thcor. Phys. 14 (1975) ^ 1, 37-54. 
[15] V. V. Dodonov and V. I. Man'ko, Invariants and correlated states of nonstationary quantum systems, in: In- 
variants and the Evolution of Nonstationary Quantum Systems, Proceedings of Lebedev Physics Institute, vol. 

183, pp. 71-181, Nauka, Moscow, 1987 [in Russian]; English translation published by Nova Science, Commack, 

New York, 1989, pp. 103-261. 
[16] P. Ehrenfest, Bemerkung iiber die angendherte Giiltigkeit der klassischen Mechanik innerhalb der Quanten- 

mechanik, Zcitschrift fiir Physik A 45 (1927), 455-457. 
[17] V. P. Ermakov, Second-order differential equations. Conditions of complete integrability , Universita Izvestia 

Kiev, Series 111 9 (1880), 1-25; see also Appl. Anal. Discrete Math. 2 (2008) #2, 123-145 for English translation 

of Ermakov's original paper. 

[18] L. D. Faddccv. Feyrman integrals for singular Lagrangians, Theoretical and Mathematical Physics 1 (1969) #1, 
3 18 [in Russian]. 

[19] P. O. Fedichcv, Yu. Kagan, G. V. Shlyapnikov and J. T. M. Walraven, Influence of nearly resonant light on the 

scattering length in low-temperature atomic gases, Phys. Rev. Lett. 77 (1996) #14, 2913-2916. 
[20] R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, New York, 1965. 
[21] S. Fliigge, Practical Quantum Mechanics, Springer- Verlag, Berlin, 1999. 

[22] T. Fujii, Sh. Matsuo, N. Hatakenaka, S. Kurihara and A. Zeilinger, Quantum circuit analog of the dynamical 

Casimir effect, Phys. Rev. B. 84 (2011) #17, 174521 (9 pages). 
[23] I. I. Gol'dman and V. D. Krivchenkov, Problems in Quantum Mechanics, Dover, New York, 1993. 
[24] J. Guerrero, and F. F. Lopez-Ruiz, V. Aldaya and F. Cossio, Harmonic states for the free particle, .1. Phys. A: 

Math. Theor. 44 (2011), 445307 (16pp); see also arXiv:1010.5525v3 [quant-ph] 1 Jul 2011. 
[25] C. H. Hagen, Scale and conformal transformations in Galilean-covariant field theory, Phys. Rev. D 5 (1972) #2, 

377-388. 

[26] G. Harari, Ya. Ben-Aryeh and Ady Mann, Propagator for the general time- dependent harmonic oscillator with 

application to an ion trap, Phys. Rev. A 84 (2011) # 6, 062104 (4 pages). 
[27] W. Heisenberg, The Physical Principles of the Quantum Theory, University of Chicago Press, Chicago, 1930; 

Dover, New York, 1949. 

[28] W. Heisenberg, Physics and Philosophy: The Revolution in Modern Science, Ruskin House, Harper and 
Row, New York, 1958 (Lectures delivered at University of St. Andrews, Scotland, Winter 1955-56); see also 
http://www.aip.org/history/heisenberg/pl3e.htm. 

[29] R. W. Henry and S. C. Glotzer, A squeezed-state pnmer. Am. J. Phys. 56 (1988) #4, 318-328. 



10 RAQUEL M. LOPEZ, SERGEI K. SUSLOV, AND JOSE M. VEGA-GUZMAN 

[30] R. Jackiw, Dynamical symmetry of the magnetic monopole, Ann. Phys. 129 (1980), 183-200. 

[31] E. G. Kalnins and W. Miller, Lie theory and separation of variables. 5. The equations iUt + Uxx = and 

iUt + U^x - c/x^U = 0, J. Math. Phys. 15 (1974) #10, 1728-1737. 
[32] Yu. S. Kivshar, T. J. Alexander and S. K. Turitsyn, Nonlinear modes of a macroscopic quantum oscillator, 

Phys. Lett. A 278 (2001) #1, 225-230. 
[33] C. Koutschan, http://hahn.la.asu.edu/~suslov/curres/index.htni; see Mathematica notebook: Koutschan.nb; 

see also http://www.risc.jku.at/people/ckoutsch/pekeris/ 
[34] C. Koutschan and D. Zeilbergcr. The 1958 Pekeris-Accad-WEIZAC Ground- Breaking Collaboration that com- 
puted Ground States of Two-Electron Atoms (and its 2010 Redux), Math. Intelligencer 33 (2011) #2, 52-57. 
[35] S. I. Kryuchkov and S. K. Suslov, On the problem of electromagnetic-field quantization, under preparation. 
[36] S. I. Kryuchkov, S. K. Suslov and J. M. Vega-Guzman, The minimum-uncertainty squeezed states for the simple 

harmonic oscillator, under preparation. 
[37] N. Lanfear and S. K. Suslov, The time- dependent Schrddinger equation, Riccati equation and Airy functions, 

arXiv:0903.3608v5 [math-ph] 22 Apr 2009. 
[38] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Nonrelativistic Theory, Pergamon Press, Oxford, 1977. 
[39] N. Lanfear, R. M. Lopez and S. K. Suslov, Exact wave functions for generalized harmonic oscillators. Journal 

of Russian Laser Research 32 (2011) #4, 352-361; sec also arXiv:11002.5119v2 [math-ph] 20 Jul 2011. 
[40] N. Lanfear, F. F. Lopez-Ruiz, S. K. Suslov and J. M. Vega-Guzman, Coherent states for generalized harmonic 

oscillators, under preparation. 
[41] P. G. L. Leach and K. Andriopoulos, The Ermakov equation: a commentary, Appl. Anal. Discrete Math. 2 

(2008) #2, 146-157. 

[42] H. R. Lewis, Jr., and W. B. Riesenfeld, An exact quantum theory of the time- dependent harmonic oscillator and 
of a charged particle in a time- dependent electrom,agnetic field, J. Math. Phys. 10 (1969) #8, 1458-1473. 

[43] R. M. Lopez and S. K. Suslov, The Cauchy problem for a forced harmonic oscillator, Revista Mexicana de 
Fisica, 55 (2009) #2, 195-215; see also arXiv:0707.1902v8 [math-ph] 27 Dec 2007. 

[44] R. M. Lopez, S. K. Suslov and J. M. Vega-Guzman, On the harmonic oscillator group, arXiv:1111.5569v2 
[math-ph] 4 Dec 2011. 

[45] R. M. Lopez, S. K. Suslov and J. M. Vega-Guzman, http://hahn.la.asu.edu/~suslov/curres/index.htm; see 

Mathematica notebook: HarmonicOscillatorGroup.nb 
[46] I. A. Malkin and V. I. Man'ko, Dynamical Symmetries and Coherent States of Quantum System, Nauka, Moscow, 
1979 [in Russian]. 

[47] E. Mcrzbachcr, Quantum Mechanics, third edition, John Wiley & Sons. Ncnv York, 1998. 

[48] W. Miller. Jr., Symmetry and Separation of Variables, Encyclopedia of Mathematics and Its Applications, Vol. 4, 

Addison- Wesley Publishing Company, Reading etc, 1977. 
[49] U. Niederer, The maximal kinematical invariance group of the free Schrddinger equations, Helv. Phys. Acta 45 

(1972), 802-810. 

[50] U. Niederer, The maximal kinematical invariance group of the harmonic oscillator, Helv. Phys. Acta 46 (1973), 
191 200. 

[51] A. F. Nikiforov, S. K. Suslov, and V. B. Uvarov, Classical Orthogonal Polynomials of a Discrete Variable, 

Springer- Verlag, Berlin, New York, 1991. 
[52] L. Pitaevskii and S. Stringari, Bose Einstein Condensation, Oxford University Press. Oxford, 2003. 
[53] S. Rosencrans, Perturbation algebra of an elliptic operator, J. Math. Anal. Appl. 56 (1976) #2, 317 329. 
[54] B. Sanborn, S. K. Suslov and L. Vinet, Dynamic invaria,nts and Berry's phase for generalized driven harmonic 

oscillators, Journal of Russian Laser Research 32 (2011) #5, 486-494; see also arXiv:1108.5144vl [math-ph] 25 

Aug 2011. 

[55] E. Schrodinger, Dei' stetige Ubergang von der Mikro-zur Makro Mechanik, Die Naturwissenshaften, 14 

(1926), 664-666; see also http://www.nobelprize.org/nobel_prizes/physics/laureates/1933/schrodinger-bio.html 

and http://www.zbp.univie.ac.at/schrodinger/euebersicht.htm. 
[56] E. Suazo and S. K. Suslov, Soliton-like solutions for nonlinear Schrddinger equation with variable quadratic 

Hamiltonians, arXiv:1010.2504v4 [math-ph] 24 Nov 2010. 
[57] E. Suazo, S. K. Suslov and J. M. Vega-Guzman, The Riccati equation and a diffusion-type equation, New York 

J. Math. 17a (2011), 225-244. 



HIDDEN SYMMETRY OF HARMONIC OSCILLATORS 



11 



[58] E. Suazo, S. K. Suslov and J. M. Vega-Guzman, The Riccati system and a diffusion-type equation, arXiv: 

1102.4630vl [math-ph] 22 Feb 2011. 
[59] S. K. Suslov, Dynamical invariants for variable quadratic Hamiltonians, Physica Scripta 81 (2010) #5, 055006 

(11 pp); see also arXiv:1002.0144v6 [math-ph] 11 Mar 2010. 
[60] S. K. Suslov, On integrability of nonautonomous nonlinear Schrddinger equations, to appear in Proc. Amer. 

Math. Soc; see also arXiv:1012.3661v3 [math-ph] 16 Apr 2011. 
[61] S. K. Suslov, The Berry phase for simple harmonic oscillators, arXiv:1112.2418v2 [quant-ph] 14 Dec 2011. 
[62] S. K. Suslov, http://hahn.la.asu.edu/~suslov/curres/index.htm; see Mathematica notebooks: BerrySummary.nb, 

Fourier. nb and Heisenberg.nb. 
[63] L. Vinet and A. Zhedanov, Representations of the Schrddinger group and matrix orthogonal polynomials , J. Phys. 

A: Math. Theor. 44 (2011) #35, 355201 (28 pages). 
[64] K. B. Wolf, On time- dependent quadratic Hamiltonians, SIAM J. Appl. Math. 40 (1981) #3, 419-431. 
[65] K-H. Yeon, K-K. Lee, Ch-I. Um, T. F. George and L. N. Pandey, Exact quantum theory of a time- dependent 

bound Hamiltonian systems, Phys. Rev. A 48 (1993) # 4, 2716-2720. 
[66] H. P. Yuen, Two-photon coherent states of the radiation field, Phys. Rev. A 13 (1976) #6, 2226-2243. 
[67] J. Q. You and F. Nori, Atomic physics and quantum optics using superconducting circuits, Nature 474 (2011), 

589-597. 

[68] A. V. Zhukov, Exact quantum theory of a time- dependent system with quadratic hamiltonian, Phys. Lett. A 256 
(1999) #5-6, 325-328. 

Mathematical, Computational and Modeling Sciences Center, Arizona State University, Tempe, 
AZ 85287-1904, U.S.A. 

E-mail address: rlopezl4@asu.edu 

School of Mathematical and Statistical Sciences & Mathematical, Computational and Modeling 
Sciences Center, Arizona State University, Tempe, AZ 85287-1804, U.S.A. 

E-mail address: sks@asu.edu 

URL: http : / /hahn . la . asu . edu/~ suslov/ index . html 

Mathematical, Computational and Modeling Sciences Center, Arizona State University, Tempe, 
AZ 85287-1904, U.S.A. 

E-mail address: jmvegaOasu.edu 



